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PREFACE 

The  work  as  reported  In  this  document  represents  the  completion 
of  the  various  tasks  as  prescribed  in  Contract  No.  DAAG29-77-G-0069, 
for  the  Army  Missile  Command,  Huntsville,  Alabama.  The  research 
is  monitored  through  the  U.S.  Army  Research  Office,  Research  Triangle 
Park,  North  Carolina  and  is  the  result  of  the  work  done  during  the 
period  Jan.  1977  - Dec.  1977.  The  program  manager  with  the  Army 
Missile  Command  was  Dr.  Donald  J.  Spring. 


SUMMARY 


The  research  effort  for  the  past  year  involved  the  development  of 
theoretical  prediction  methods  for  the  aerodynamic  loading  on  a wing  with 
a full  span  elevon.  The  methods  are  based  on  lifting  surface  kernel 
function  formulations  in  both  subsonic  and  supersonic  potential  flow. 

The  unique  idea  in  both  cases  is  the  closed  form-finite  summation  manner 
in  which  the  kernel  function  integral  is  solved.  This  method  of  solution 
avoids  Mangier-  and  Cauchy-type  singularity  problems,  encountered  in 
classical  numerical  integration  approaches,  and  leads  to  stable,  rapidly 
convergent  solutions.  In  subsonic  flow,  an  existing  kernel  function 
method  for  planar  wings  was  modified  by  adding  an  assumed  pressure  loading 
function  to  account  for  the  presence  of  the  elevon.  The  assumed  pressure 
loading  distribution  led  to  exact  closed  form  solutions  for  section  and 
total  coefficients  on  the  wing;  however,  for  the  case  of  the  deflected 
elevon,  some  numerical  integration  procedures  were  required.  Results  of 
these  computations  agree  very  well  with  experimental  data.  In  the  super- 
sonic Mach  number  regime,  a lifting  surface  kernel  function  method  similar 
to  the  subsonic  approach  was  developed  for  the  planar  wing  case,  but  with 
appropriate  Mach  cone  regions  of  integration  taken  into  account.  Various 
assumed  pressure  loading  functions,  all  weighted  by  exact  theoretical 
results,  were  required  for  different  wing  shapes.  Numerical  results 
produced  stable  (nonosclllatory)  solutions  which  agreed  well  with 
experimental  data,  even  for  very  low  aspect  ratio  triangular  wings. 
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INTRODUCTION 
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I The  problem  of  predicting  the  wing  and  control  surface  loading  and 

associated  hinge  moments  of  a missile  in  both  subsonic  and  supersonic 

I 

flight  has  been  the  subject  of  the  present  research  for  the  last  year 
(since  January  1977).  More  specifically,  the  loading  for  a "deformed  wing" 

I 

I 

has  been  studied  and  preliminary  results  have  been  obtained.  Problems  of 
this  nature  have  been  extensively  studied  (Refs.  1-35)  and  can  be  approached 
in  three  fundamental  ways.  The  first  and  perhaps  the  most  expensive  is  the 
experimental  approach.  Wind  tunnel  tests  have  been  used  extensively  for 
problems  of  this  nature  for  many  years.  For  obvious  reasons,  this  procedure 
is  rather  expensive  and  long  lead  times  are  required.  Additionally,  exten- 
sion of  the  data  to  other  configurations  is  somewhat  questionable  and  can 
sometimes  lead  to  erroneous  conclusions. 

The  second  approach  is  a finite  element  approach.  The  flow  field  is 
subdivided  into  small  finite  elements  and  the  field  properties  are  com- 
puted at  each  element.  For  most  real  problems,  the  element  size  must  be 
significantly  smaller  than  some  characteristic  dimension  of  the  flow  field 
and  consequently  elements  numbering  in  the  thousands  and  tens  of  thousands 
are  generally  required  even  for  two-dimensional  analyses.  For  three- 
dimensional  problems,  another  order  of  magnitude  must  be  added  which  would 
require  very  long  computer  run  times.  The  advantages  of  such  approaches 
are  good  accuracy  and  usually  good  detailed  information  concerning  the 
flow  field,  aerodynamic  loads  and  field  properties.  The  scheme  is  not 
limited  to  inviscid  flow  nor  to  small  perturbations.  However,  it  was  felt 


1 


2 


that  the  disadvantages  of  long  computer  run  times,  high  cost,  and  re- 
striction for  full  three-dimensional  flows  prohibit  the  use  of  finite 
element  or  finite  difference  approaches.  .• 

The  third  approach  is  fundamentally  founded  in  linearized  potential 
flow  theory  and  is  the  one  taken  in  the  present  research  work.  Many 
techniques  for  solving  the  potential  flow  problem  are  available  and  two 
are  basically  used  in  wing  loading  computations.  The  first  is  based  on 
vortex  lattice  theory  and  is  used  extensively  for  problems  concerning 
wings  and  other  lifting  bodies.  However,  it  too  requires  large  computer 
storage  (normally)  and  relatively  long  run  times.  Additionally,  vortex 
lattice  techniques  fail  to  accurately  predict  the  wing  loading  near  trail- 
ing edges  (Ref.  37)  and  tend  to  produce  oscillatory  pressure  distributions 
in  supersonic  flow  (Ref.  31). 

For  the  present  research,  the  lifting  surface  (kernel  function) 
formulation  is  used.  This  approach  is  known  to  give  good  agreement  between 
theory  and  experiment  but  past  formulations  have  required  integration  over 
a second  order  singularity.  This  problem  has,  however,  been  overcome  in 
the  present  work  by  removing  the  pressure  loading  term  in  the  equations 
which  allows  direct  closed  form  integration  of  the  remaining  terms. 

The  coding  described  above  has  been  used  to  determine  the  loading  on 
wings  and  elevons  and  hinge  moment  calculations  have  been  made  for  these 
configurations.  The  analysis  is  also  applicable  to  the  "detached"  eleven 
case  in  which  a gap*  exists  between  the  eleven  leading  edge  and  the  wing 
trailing  edge. 

*Under  the  assumption  that  the  gap  is  less  than  local  wing  chord  and 
wing  wake  has  no  significant  rollup. 
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The  use  of  the  kernel  function  formulation  as  described  in  this  report 
retains  the  advantages  of  the  theory  while  eliminating  or  minimizing 
some  of  its  inherent  problems.  Lifting  line,  thin  airfoil  and  slender  wing 
theory  results  are  used  to  develop  general  functional  forms  for  the 
distributions  and  the  basic  integral  equation  is  expressed  in  terms  of 
vorticity  which  satisfies  the  linearized  compressible  potential  equation. 
The  functional  forms  for  the  spanwise  and  chordwise  loadings  are  by  no 
means  based  on  empirical  results  from  experimental  data  but  are  logically 
deduced  from  exact  theories  as  mentioned  above.  In  fact,  the  thin  air- 
foil and  lifting  line  theories  can  be  shown  to  be  special  limiting  cases 
of  the  general  lifting  surface  integral  (Ref.  34).  Evaluation  of  the 
integral  results  are  very  similar  to  the  vortex  lattice  equations.  Use  of 
the  Multhopp  spacing  (for  subsonic  flow)  for  control  points  in  subsonic 
flow,  in  conjunction  with  the  above,  produces  accurate  flow  tangency  over 
the  entire  wing  without  the  use  of  least-squares  formulation. 

In  order  to  provide  a complete  formulation  of  the  problem  (in  one 
document)  the  basic  elements  of  the  entire  derivation  of  the  potential 
equation  is  presented  in  the  Appendix. 

For  supersonic  flow,  the  basic  approach  uses  potential  flow  theory 
as  in  the  subsonic  case.  This  allows  the  general  problem  to  be  modeled  as 
a superposition  of  the  various  singularities  with  the  appropriate  boundary 
conditions.  There  is  an  additional  requirement  that  the  flow  remain 
attached  at  every  point  on  the  lifting  surface.  The  pressure  distribution 
on  the  wing  surface  is  found  by  solving  the  supersonic  kernel  function, 
which  relates  the  pressure  distribution  over  the  surface  to  the  local 
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surface  velocities.  The  details  of  the  theory  and  numerical  method  will 
be  given  in  the  Supersonic  Analysis  section  following  a brief  review  of 
various  current  methods  used  in  supersonic  potential  flow  problems. 


SUBSONIC  ANALYSIS 

From  the  Appendix  it  is  shown  that  the  potential  equation  for  a 
planar  lifting  surface  in  compressible  flow  is 


(|i(x 


= h If 


(y-y^)^  + z' 


Z 1 + 


(x-x^) 


V (x-x  (y-y  )2+b2z2 


I dx 

J 


dy. 


Using  the  relation  between  the  surface  vorticity  strength  and  the 
pressure  loading  coefficient  y = 40^/2  and  employing  the  dimensionless 
spanwise  and  chordwise  variables  n and  the  nondimens ional  normalwash 
in  the  plane  of  the  sheet  due  to  the  sheet  itself  is 


= i ff  i -vd 


1 + 


(x-x^) 


(y-yo)" 


Vcx-x 


)2+g2(y_y 

O O 


r-  1 dy^ 

)2  J 


(2) 


where 

n * y/ (b/2) 

and 

C = (x-x^g)/c(n) 

If  AC  is  expressed  as  the  proper  combination  of  spanwise  and  chord- 
P 

wise  variables,  the  behavior  along  a chord  or  across  the  span  can  be 
examined  Independently.  Since  the  loading  along  each  chord  should  behave 


5 


in  the  classical  two-dimensional  manner  at  the  leading  and  trailing 
edges,  thin  airfoil  theory  results  can  be  used  to  define  AC  as 


AC  (C.n) 
P 


A (n)c" 
n 


(3) 


The  first  term  (n=0)  produces  the  required  leading  edge  singularity  (at 
C=0)  and  each  term  satisfies  the  Kutta  condition  at  the  trailing  edge. 

Lifting  line  theory  indicates  that  the  spanwise  coefficients  should 
be  elliptic  in  nature  for  finite  values  of  taper  ratio,  while  for  zero 
taper  ratio  the  required  behavior  is  that  of  slender  wing  theory.  Intro- 
ducing the  nondimensional  spanwise  variable,  6,  where  n = cos  0,  a 
functional  form  satisfying  both  of  these  requirements  is 


A (n) 
n 


1 

c(n) 


M 

I 


m=0 


B sin(2m+l)9 
nm 


(4) 


The  B 's  are  constant  coefficients  and  the  c(n)  term  is  the  local  chord, 
nm 

To  check  the  validity  of  Equation  (4),  consider  a delta  wing.  For  m=0, 
we  then  have 


An(n)  = 


cj^(l-n) 


[B^j^  sin0] 


However,  sin  0 may  be  written  as,  sin  0 

A (n)  - 


1-n^ , and  hence 


which  is  the  required  behavior  for  slender  wings.  For  the  case  of  a 
straight  nontapered  wing,  the  m-O  term  reduces  to 


6 


A„(n)  = 

n c 


nl 


which  is  the  required  elliptic  spanwise  loading  term. 

Other  terms  in  Equation  (4)  corresponding  to  m = 1,2,3,...,  merely 
form  corrections  of  a Fourier  Series  nature  to  the  basic  elliptic  term 
which  dominates  the  solution.  It  may  be  argued  that  other  functional 
forms  for  the  spanwise  loading  could  be  used  and  that  the  one  presented 
in  Equation  (4)  is  not  unique.  Indeed  this  is  true  and  other  forms  have 
been  used  but  the  trigometric  functional  sin  (2m+l)0  seems  to  produce 
the  better  results  at  least  for  subsonic  flow.  However,  any  solution  which 
is  used  must  be  dominated  by  the  basic  elliptic  loading  terms. 

If  Equations  (3)  and  (4)  are  now  combined,  the  complete  expression 
for  the  pressure  loading  coefficient  becomes 


M 

I 

m=0 


(5) 


Spanwise 


Chordwise 


Returning  to  the  evaluation  of  the  downwash,  note  that  over  a suf- 
ficiently small  element  of  the  planform  surface  AC^  is  approximately 
constant,  and  so  can  be  taken  outside  the  integral.  This  procedure  is 
analogous  to  the  constant  panel  loading  assumption  which  is  fundamental 
to  vortex  lattice  theory.  The  downwash  due  to  this  small  element  then 
becomes 


Aw(x,y) 


ACp(C,n) 

8it 


22  / (x-x„) 

f / — ^ (if  - 

i.  ' V(X-X  )2+b2  )2 

1 r o o 


(6) 
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where  AC^(5,n)  is  evaluated  at  the  centroid  of  the  element.  The  integral 
thus  obtained  can  be  evaluated  in  closed  form  with  the  result 


Aw(x,y) 


AC  a,n) 


[K(x2,y2)  -^(x^.yj^)  -KCx^.y^)  +K.(x^,y^)] 


where 


(x-x  ) +W(x-x  )2  + 02  ( )2 

K(x  ,y  ) = 9 1 2 1 ^.o. 

° ° f \ 

(y-yo) 


2 j.  q2/'„_„  '»2i  (8) 


+ B logg[B(y-yQ)  +'y(x-x^)^  + B^Cy-y^)^] 


If  the  entire  planform  is  divided  into  similar  small  elements,  the 
total  downwash  at  any  point  on  the  wing  surface  may  be  written  as 

w(x,y)  = [ Aw(x,y) 

S 

where  Equations  (5)  and  (7)  are  used  to  evaluate  Aw. 

The  nondlmensional  tangency  condition  for  small  angles  of  attack  is 


w(x,y)  + sin  a(x,y)  = 0 


which  defines  the  load  distribution  in  terms  of  the  physical  requirements 
of  no  flow  through  the  wing.  Choosing  spanwise  and  chordwise  control 
points  according  to  a cosine  distribution 


^ -^2 


(1  - cos  ^),  i = 1,2,...N+1 


“ I ^ ■ 1.2.... M+1 
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and  evaluating  Equation  (10)  at  each  control  point  gives: 


J Aw(x  ,y.)  + sin  a(x.,y,)  = 0 
c ^ J 1 J 


Equation  (12),  together  with  the  defining  relations  (5)  and  (7),  con- 
stitutes a set  of  simultaneous  linear  algebraic  equations  for  NxM 

unknown  loading  coefficients  B . Solving  the  system  for  these  coef- 

nm 

ficlents,  the  pressure  loading  at  any  point  on  the  wing  is  then  available 
from  Equation  (5). 

The  spanwise  section  loading  is  given  by 


X 

cc^(yo)  = c(y^) J ACp(C.n^)dC 


This  integration  can  be  done  in  closed  form.  Defining 


.n  = 0,l,...,N 


and  integrating  yields 


I = ^ 
o 2 


T - 2n-l 

^n  2n+2  ^n-1  ’ " “ 


For  convenience,  let 
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then  after  substitution  into  Equation  (13)  from  Equations  (7) , (15)  and 
(16),  successively  yields 


M 


m=0 


Similar  manipulations  lead  to 


(17) 


cc  (y  ) M rsin(2m+l)e 

^^d^yo>  = — I (2m+l)E^  — °l 

do  8 m sin  9^ 
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[sin(2m+l)6  I 

®o  J 


(18) 


and 


M N 

I I B I sin(2m+l)e 
mo  o'  ^ nm  n+1  o 

m=0  n=0 


(19) 


for  the  section  drag  and  pitching  moment,  respectively.  The  total  lift 
and  induced  drag  are  found  by  evaluating  the  standard  integral  forms 


b/2 

“ s f 


(20) 
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and 


b/2 

S “ I / "‘^d^yo^'^yc 
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(21) 


Using  the  spanwlse  angular  variable  0 and  Equations  (17)  and  (18)  in  the 


above  results  in 
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ttE 


C,  = 


and 


L 2S 
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(22) 
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(23) 


which  are  the  classical  forms  similar  to  those  obtained  from  Glauert's 
lifting  line  analysis. 

The  total  pitching  moment  about  the  wing  root  chord  leading  edge  is 
found  by 


b/2 

ir  / 


- J ccj^(y)x(y)dy 
-b/2 


(24) 


as 


Referring  to  Figure  1,  it  is  seen  that  Equation  (24)  may  be  written 
1 

*^m  - / (y  tan  A + x )dy 

J ^ CD  ■' 


or 


C = - i- 
m 

Sc 


r (cc.  y tan  A + cc  )dv 

“le 

0 


(25) 


Using  Equations  (16),  (17)  and  (18),  Equation  (25)  results  in 


C =. 


Sc 


- ^m  " 2 ‘^r 


m=0 


M 

tan  A I E G 1 

r\  ® m 

in“0 


(26) 
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N 7r/2 

where  F = T B I and  G = / sin  26  sin  (2in+l)0  d0. 

n=0  0 

Note  that  the  first  term  in  Equation  (26)  is  the  moment  contribution 
due  to  wing  taper,  and  the  third  term  is  the  contribution  due  to  leading 
edge  sweep.  The  second  term  looks  very  much  like  the  classical  Glauert 
term  from  lifting  line  theory. 

The  preceeding  equations  constitute  a set  which  determines  the 
fundamental  parameters  of  interest;  lift,  induced  drag,  and  pitching 
moment  for  a thin  wing  in  subsonic  compressible  flow.  Before  proceeding 
with  extensions  of  the  analysis  to  elevons,  comments  on  the  accuracy  of 
the  theory  are  in  order. 

The  success  or  failure  of  any  theory  rests  on  its  agreement  with 
experiment;  however,  before  comparing  with  experiment  the  accuracy  of 
the  flowfield  model  should  first  be  determined.  The  tangency  of  the  flow 
to  the  wing  surface  indicates  both  how  well  the  AC^  distribution  has  been 
represented  and  to  what  degree  the  wing  has  been  replaced  with  a solid 
surface.  Obviously,  low  aspect  ratio  delta  wings  are  the  severest  test 
for  planar  lifting  surface  theories  due  to  the  decreasing  area  and  large 
ACp  gradients  near  the  wing  tips.  Figure  2 shows  the  distribution  of  a-v, 
or  net  normalwash,  over  the  surface  of  a delta  wing  of  aspect  ratio  2. 

From  the  present  analysis,  less  than  two  percent  deviation  (0.02a)  from  the 
desired  zero  value  is  obtained  except  near  the  leading  edge,  which  indicates 
that  most  of  the  wing  has  Indeed  been  replaced  by  a solid  surface.  The 
behavior  near  the  leading  edge  is  due  to  computer  difficulties  in 
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evaluating  / (1-C)/C  near  i=0.  Compared  to  the  flowfields  in  vortex 

lattice  theory,  even  the  errors  near  the  leading  edge  and  the  tips  would 

be  considered  minor. 

i AC  Distribution 

' — P 

In  addition  to  the  net  normalwash  distribution  shown  in  Figure  2,  a 
comparison  was  made  with  experimental  data  for  the  pressure  distribution 
over  the  delta  wing  (Wick  wing)  sketched  in  Figure  3.  Results  of  one  of 
these  comparisons  is  shown  in  Figure  4.  The  experimental  data  extracted 
from  Ref.  (38)  agrees  very  well  with  the  theory  except  for  the  differences 
near  the  leading  edge.  These  deviations  are  easily  accounted  for  due  to 
the  wing  thickness. 

Lift  and  Pitching  Moment  Results 

To  illustrate  the  utility  and  accuracy  of  the  lift  and  pitching  moment 
calculations  of  the  analysis,  several  wing  geometries  were  evaluated.  In 
Figure  5,  two  wings  with  the  same  leading  edge  sweep  angle,  40°,  but 
different  taper  ratios  were  evaluated.  The  spanwlse  lift  distribution 
agrees  very  well  with  experimental  data  and  the  spanwlse  center  of  pressure 
and  hence  the  section  pitching  moment  is  essentially  within  experimental 
error. 

The  measured  and  predicted  C^^  's  for  low  aspect  ratio  delta  wings  are 

a 

shown  in  Figure  6.  The  theoretical  curves  properly  approach  the  slender 
wing  theory  results  for  very  small  aspect  ratios,  and  tend  toward 
2t[  as  AR 


Figure  7 shows  the  effect  of  Mach  number  on  the  theoretical  and 
experimental  lift  curve  slopes  of  two  thin  wings.  The  agreement  between 
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Figure  A.  Planfonn  Pressure  Distribution 
on  Delta  Wing  (M  “ -13) 
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Aspect  Ratio 
(delta  wings) 


Experiment  ^ O 

Present  Method 

Theory  Numerical  Integration 

— Slender  Wing 

Figure  6.  Lift  Curve  Slopes  for  Delta  Wings. 
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theory  and  experiment  is  excellent  for  Mach  numbers  less  than  approximately 
0.85.  The  good  agreement  above  M = 0.6  is  due  to  several  factors: 

(1)  both  wings  are  ideally  suited  to  a planar  wing  analysis  since  the 
thickness-to-chord  ratios  of  both  wing  sections  are  less  than  five  percent; 

(2)  the  experimental  C 's  were  evaluated  at  low  angles  of  attack;  (3)  and, 

a 

most  significantly,  the  compressibility  factor  B enters  the  theoretical 
equations  in  the  same  manner  as  in  slender  body  theory  (Ref.  3);  therefore, 
based  on  slender  body  theory  results,  the  true  transonic  region  should  be 
much  smaller  than  that  usually  assumed  in  wing  theory. 

Thickness  Effects 

Although  thickness  is  not  considered  to  be  a major  analytical  objective 
in  the  present  analysis,  thickness  effects  were  considered.  As  an  example, 
the  wing  sketched  in  Figure  3 is  relatively  thick  at  least  along  the  center- 
line.  Section  lift  coefficients  for  this  wing  at  various  angles  of  attack 
are  presented  in  Figure  8.  In  this  figure,  no  corrections  have  been  made 
for  thickness  and  still  the  agreement  between  theory  and  experiment  is 
relatively  good.  The  corresponding  center  of  pressure  (and  hence  pitching 
moment)  calculation  is  shown  in  Figure  9 with  excellent  agreement  with 
experimental  data. 

To  account  for  thickness,  two  approaches  were  considered.  First,  the 
existing  lifting  surface  program  was  modified  so  that  boundary  conditions 
were  met  on  the  wing  surface  as  opposed  to  the  "thin  wing"  planar  surface. 
Results  of  this  analysis  is  shown  in  Figure  10. 

In  the  second  approach,  the  wing  subpanels  were  overlaid  with  source 
panels.  This  approach  proved  to  be  rather  lengthy  and  complex  but  results 
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Figure  10.  Spanwlse  Loading  for  Thick  Wing  (Wick) 
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were  obtained  which  indicated  that  this  was  indeed  a feasible  approach. 
However,  because  of  its  long  computer  run  time  and  large  storage,  the 
procedure  was  not  Included  in  the  main  analysis  and  computer  program. 

Elevon  Loading  and  Hinge  Moments 

When  a control  surface  (elevon)  is  added  to  a wing  at  the  trailing 
edge,  obviously  the  loading  distributions  on  the  wing  will  change.  For 
a deflected  elevon,  a discontinuity  exists  in  the  wing  surface  slope  and 
theoretically  the  velocity  and  hence  the  pressure  differential  is  Infinite 
at  the  hinge  line.  For  the  present  anlaysis  it  is  assumed  that  the 
hinge  line  for  the  elevon  is  at  the  elevon  leading  edge  or  wing  trailing 
edge.  Additionally,  it  is  assumed  that  the  elevon  span  is  that  of  the 
wing  and  that  the  hinge  line  is  perpendicular  to  the  wing  centerline  (root 
chord) . The  analysis  under  these  assumptions  parallels  that  of  the  basic 
wing  with  the  exception  of  the  loading  functions. 

Two  approaches  were  considered  in  the  analysis  of  the  wing  loading 
with  a deflected  elevon.  Glauert  (39)  used  a single  function  Fourier 
series  to  represent  the  two-dimensional  chordwise  load  distribution.  How- 
ever, a very  large  number  of  terms  are  required  in  order  to  adequately 
represent  the  discontinuity  at  the  elevon  leading  edge.  Since  the  tech- 
niques used  in  this  analysis  is  somewhat  "self-correcting"  a simpler  single 
function  was  initially  used  to  represent  the  wing  loading  as 


AC 

P 


F(n)c" 


(27) 
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Results  using  a function  of  this  type.  Equation  (27),  proved  to  be 
very  successful  in  predicting  the  AC^  distribution;  however,  integration 
of  the  pressure  distribution  over  the  wing  is  rather  difficult  and  must 
be  done  numerically.  Schematically,  the  loading  for  the  wing  using 
Equation  (27)  is  shown  in  Figure  11a.  Because  of  the  difficulty  in  inte- 
gration (no  closed  formed  solution)  the  single  function  approach  was 
temporarily  shelved  and  a second  method  was  tried. 

In  the  second  approach,  two  functions  were  used;  one  to  represent 
the  wing  loading  and  one  to  represent  the  eleven  loading.  For  the  wing 
loading  term,  the  leading  edge  singularity  is  modeled  as  before  but  the 
Kutta  condition  is  not  satisfied  at  the  trailing  edge  of  the  wing.  The 
eleven  loading  term  has  the  required  leading  edge  singularity  and  the 
Kutta  condition  is  satisfied  at  the  trailing  edge  of  the  elevon.  The 
loading  function  for  the  two  functions'  approach  is  shown  in  Figure  11b. 

The  vorticity  (or  AC^)  at  any  point  on  the  wing  or  elevon  is  de- 
termined by  solving  the  set  of  equations  as  defined  by  Equation  (12); 
however,  loading  of  the  matrix  is  somewhat  more  complex.  For  points  on 
the  wing,  the  boundary  condition  is  (as  before  for  the  wing  above)  sin  a. 
However,  on  the  elevon,  the  boundary  condition  is  sin  (a+6)  where  6 is 
the  elevon  deflection  angle  relative  to  the  wing. 

In  order  to  clearly  illustrate  the  method  by  which  the  unknown  co- 
efficients, B , are  determined,  consider  an  example  wing  as  shown  in 
nm 

Figure  12a.  Suppose  4 control  point  locations  are  placed  on  the  wing 
in  the  chordwise  direction,  A on  the  elevon  in  the  spanwlse  direction  and 
2 on  the  elevon  in  the  chordwise  direction.  Schematically,  the  control 
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Figure  12.  Typical  Control  Point  Locations  and  Matrix  Schem.tlc 
for  Wlng-Elevon  Configuration. 
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The  matrix  for  determining  the  unknown  B's  is  shown  in  Figure  12b. 

Note  that  any  one  row  of  the  coefficient  matrix  represents  the  velocities 

Induced  at  one  control  point.  The  numbering  scheme  for  the  control  point 

locations  is  also  shown  in  Figure  12,  The  B 's  are  the  coefficients 

nm 

associated  with  the  wing  loading  function  and  the  ^pq'®  associated 
with  the  loading  function  for  the  elevon.  Once  the  B's  have  been  de- 
termined, the  loading  on  the  wing  or  elevon  may  be  computed  independently 
as  shown  in  the  next  section. 

Sectional  Loads  and  Moments 

The  sectional  lift  for  a wing  with  a deflected  elevon  is 
x=x,.  X=r 


cc 


r “ 

,(y  ) = / AC  (C,n)  d + / Ac_  (c,n)dx 

O J p^  X y Pg 

x=x^ 


(28) 
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where  the  total  chord  length  is 


‘^To"  ‘^E 


As  before  (for  the  wing  alone)  define 


(29) 


C = x/c, 


To 


(30) 


Then  the  first  Integral  for  the  wing  becomes 
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For  the  second  Integral  define 


(32) 
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and  Equation  (28)  may  then  be  written  as 


cc  = c 
I Xo 


J AC  d C + Cg  r AC  (Cj,.n)dC 
=0 


-n  /1-C 
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Now,  for  the  wing 

Nw  M 

= 'l  , B sin(2m+l)e  C 

Pw  n=0  m^O  " V ^ 

In  order  to  evaluate  the  first  integral  in  Equation  33,  define 
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Integration  of  Equation  (35)  yields 


where 
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For  other  values  of  n.  Equation  (35)  becomes 


^wn  = 'n-1 
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For  the  second  Integral  in  Equation  (33),  the  results  are  for  the  eleven 
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and 
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Hence, 
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The  total  lift  for  the  wing  then  becomes 
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cc^(n)  d n = l/  cc^  dn 
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Evaluation  of  Equation  (41)  is  not  easily  done  analytically  and 

hence  the  integral  is  done  numerically  using  the  trapezoidal  rule.  The 

reason  for  the  difficulty  lies  in  the  fact  that  for  a general  wing  plan- 

form  shape,  C is  a function  of  n and  hence  even  for  n=0  the  integral 
n 

is  very  complex.  For  the  elevon  alone,  one  could  in  fact  find  the  closed 
form  solution  for  the  total  lift  (on  the  elevon) . 

The  pitching  moment  for  the  entire  wing  and  elevon  is  found  by 
evaluating 


1 


where  cc  is  the  section  pitching  moment  about  the  y axis  as  shown  in 
mrc  r o 

Figure  1.  Equation  (42)  is  also  evaluated  numerically. 


31 


Of  particular  interest  where  the  eleven  is  concerned  is  the  compu- 
tation of  the  eleven  hinge  moment.  This  could  be  extracted  numerically 
from  Equation  (42);  however,  it  may  be  computed  in  closed  form.  The 
moment  coefficient  about  the  hinge  line  is  (See  Fig,  13) 


~ / cc  (y)  X.  (y)dy 

-b/2  E ® 


where 


The  ACp  distribution  for  the  eleven  is  (see  Equation  (34)  for  the  wing), 
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Following  the  general  procedure  as  outlined  for  the  wing,  the  eleven  hinge 
moment  is  evaluated  by  substituting  Equation  (46)  into  Equation  (44)  and 


the  result  is 
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where 


and 
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Finally,  Equation  (47)  is  substituted  into  Equation  (44)  and  the  result 
is 
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where  G = / sin  26  sin  [(2q+l)e]d9 
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Results  of  Elevon  Loading  Computations 

At  the  present  time  the  elevon-wing  analysis  and  results  are  in 
the  beginning  stages  and  very  little  experimental  data  seems  to  be  available 
for  comparison  purposes.  However,  "some"  experimental  data  was  found  for 
four  wings  ranging  in  aspect  ratios  from  2 to  6 and  results  are  shown  in 
Figure  14. 

In  Figure  14,  experimental  values  of  C and  C are  compared  with  the 

a ^6 

present  theory  with  satisfactory  results.  The  theory  utilizes  the  two- 
function  approach  and  is  within  experimental  error  of  the  data. 
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In  sununary,  the  objectives  for  the  research  effort  for  the  first 
year  have  been  accomplished.  Excellent  results  have  been  obtained  in 
some  areas  and  areas  where  further  concentrated  efforts  are  required  have 
been  Identified. 


SUPERSONIC  ANALYSIS 

Current  Methods 

There  are  several  methods  available  for  solving  supersonic  potential 
flow  problems.  These  include  exact  theories,  discrete  element  methods, 
kernel  function  methods,  and  numerical  approached  (finite  element  or 
finite  difference)  schemes.  All  potential  flow  methods  may  be  grouped  in 
one  of  two  categories;  linear  theory  or  full  potential  (nonlinear  numer- 
ical approach) . 

Numerical  methods  for  solving  the  nonlinear  potential  equations,  such 
as  finite  difference  techniques,  provide  the  most  accurate  solutions  avail- 
able. These  methods  were  not  considered  for  the  current  problem,  however, 
due  to  the  high  cost  and  excessive  computation  times  usually  required. 

Linearized  theory  methods,  based  on  the  superposition  of  elemental 
solutions  of  the  linearized  differential  equation  governing  the  flow, 
have  been  employed  extensively  in  analysis  of  supersonic  aircraft.  Exact 
theories,  discrete  element  methods  and  kernel  function  methods  are  all  in 
the  realm  of  linear  theory. 

Several  exact  linear  theory  methods  are  available,  such  as  Slender 
Body  Theory  (44),  Conical  Flow  Theory  (43)  and  Eward's  method  (28,45). 
Slender  Body  Theory  uses  two-dimensional  crossflow  and  slender  wing-body 
ideas  to  generate  general  forces  and  moments  for  an  entire  configuration. 
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Details  of  the  pressure  loading  are  not  available  and  general  deformed 
wings  cannot  be  handled.  Conical  flow  theory  gives  exact  equations  for 
over  thin  flat  wings  with  subsonic  leading  edges  at  constant  angle 
of  attack.  Eward's  theory  involves  source  sheet  Integrals  representing 
wings  with  supersonic  leading  and  trailing  edges.  Exact  expressions  for 
ACp  based  on  Eward's  theory  have  been  derived  for  constant  angle  of 
attack,  constant  roll  rate,  and  constant  pitch  rate.  The  general  formu- 
lation is  also  the  basis  for  several  discrete  element  methods.  Unfortun- 
ately, none  of  these  methods  give  results  which  are  generally  applicable 
to  a deformed  wing  study.  The  solutions  do,  however,  give  insight  into 
flow  patterns  and  load  distribution  trends,  and  can  be  used  as  limiting 
check  cases  for  a deformed  wing  method. 

Discrete  element  methods  are  computerized  numerical  techniques  which 
use  elemental  source  and  vortex  solutions,  in  conjunction  with  the  super- 
position principle,  to  model  the  flow  over  general  configurations.  In 
these  methods  the  wing  is  subdivided  into  a mesh  or  grid,  each  panel  of 
which  is  a source  or  vortex  sheet  of  constant,  but  usually  unknown, 
strength.  The  unknowns  are  found  either  by  inverting  a large  influence 
coefficient  matrix,  representing  the  flow  tangency  condition  at  selected 
points  of  the  wing,  or  by  a downstream  marching  technique.  Current  methods 
are  available  which  can  solve  problems  with  thousands  of  panels;  however, 
the  cost  per  solution  increases  with  increasing  number  of  panels.  Examples 
of  the  discrete  element  method  include  the  approaches  of  Etkin  (46)  and 
Purvis  (47),  which  are  based  on  Eward's  Theory,  the  well  known  Woodward 
technique  (7),  which  uses  source  and  vortex  elements,  the  "Mach  box"  method 
of  Pines  and  Dugundjl  (48),  and  the  panel  method  of  Appa  and  Smith  (49). 
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Supersonic  kernel  function  methods  have  also  been  developed  with 
varying  degrees  of  success.  Some  analysts,  such  as  Cunningham  (32)  and 
Curtis  and  Llngard  (50) , use  assumed  functions  to  solve  the  wing  problem 
in  the  same  manner  as  that  employed  in  subsonic  kernel  function  methods. 
Others,  such  as  Miranda  (51)  and  Carlson,  et  al,  (52,31,53),  solve  the 
kernel  Integral  in  a manner  similar  to  the  discrete  element  approaches. 

The  most  serious  disadvantage  of  the  kernel  function  method  is  that  plan- 
forms  with  discontinuous  pressure  distributions  are  very  difficult  to 
treat.  Since  the  assumed  pressure  functions  are  smooth,  the  convergence 
to  a discontinuous  distribution  is  very  slow.  When  discrete  elements  are 
used,  high  frequency  disturbances  propagate  along  Mach  lines  and  may 
increase  in  amplitude.  Cunningham  (32)  avoided  the  former  problem  by 
using  exact  theory  weighting  functions,  which  contain  the  appropriate 
discontinuities.  In  the  latter  case,  Carlson,  et  al  (31,53,28),  employed 
a "terminal  smoothing"  procedure  to  eliminate  the  undesired  wiggles. 

The  present  method  proposes  to  combine  the  assumed  pressure  function 
approach  for  solving  the  supersonic  kernel  integral  with  the  discrete 
element  technique.  Results  will  show  that  this  approach  retains  the  speed 
and  accuracy  advantages  of  the  kernel  function  method  without  encountering 
any  of  the  previously  mentioned  problems. 

Theoretical  Approach 

A thin  planar  lifting  surface  in  supersonic  flow  may  be  represented 


mathematically  by  a finite  sheet  of  elemental  horseshoe  vortices  in 

the  same  manner  as  in  subsonic  flow  (Ref.  21).  The  nondimensional  velocity 

potential  for  such  a surface  is  (See  Appendix) 
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V (y-y„)2+z2  L 


(x-Xo) 


Vcx-X  )2-62(y-y  )2_^2^2 


]dx  dy 
o 


(51) 


The  variables  of  integration  x^,y^  are  defined  only  over  the  planform 
surface,  S,  and  *^he  unknown  vorticity.  Equation  (51)  differs 

from  the  subsonic  result  by  (1)  the  term  in  brackets,  (2)  the  factor  Zir 
Instead  of  4tt,  and  (3)  the  region  of  integration  includes  only  that 
portion  of  the  wing  surface  S which  is  contained  in  the  Mach  forecone  from 
point  (x,y,z). 

The  relation  between  the  pressure  loading  coefficient  AC  and  the 

P 

vorticity  is  the  same  as  that  for  subsonic  flow.  Applying  the  Kutta- 
Joukowskl  Theorem  at  a point  on  the  sheet  gives 


AC  (x  ,y  ) 


(52) 


where  AC  = C - C 

^ ^lower  '^upper 

The  nondimenslonal  perturbation  velocity  w,  or  downwash,  which  the 
lifting  surface  induces  normal  to  itself  is  found  by  differentiating 
Equation  (51)  with  respect  to  z and  taking  the  limit  as  z approaches 
zero.  This  process  results  in 


w(x,y) 
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AC  (x  ,y  ) 
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(x-x  ) 
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dx  dy 
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(53) 


which  is  the  well  known  supersonic  kernel  function  integral.  The  down- 
wash  that  a lifting  surface  in  supersonic  flow  induces  on  Itself  is 
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unique  in  that  the  downwash  point  is  located  on  the  boundary  of  the 
region  of  Integration.  The  integral  gives  the  appearance  of  being 
improper  due  to  the  singularity  at  y=y^  within  the  region  of  integration. 
This  Integrand  is,  as  previously  shown,  the  limiting  form  in  the  z=0  plane 
of  a more  general  non-singular  integrand.  Consequently,  the  Integral  would 
be  expected  to  yield  a finite  result  for  w just  as  in  the  subsonic  case. 
There  also  appears  to  be  an  upwash  singularity  along  Mach  waves  through 
x^,y^  due  to  the  square  root  term  in  the  denominator.  As  noted  in 
Reference  31:  "...this  same  upwash  field,  with  the  large  values  of  upwash 
near  the  Mach  cone  limits,  is  also  largely  responsible  for  the  difficulties 
in  representing  supersonic-flow  phenomena  by  means  of  finite-element 
techniques."  This  problem  does  not  appear  in  the  current  method,  which 
contains  finite-element  type  terms,  and  in  fact,  the  influence  of  a 
discrete  panel  along  its  Mach  waves  will  be  shown  to  be  a finite  value. 

Equation  (53)  is  the  starting  point  adopted  by  most  authors  (31,32) 
for  the  planar  wing  case. 

The  first  step  in  solving  the  Integral  equation  is  to  determine  the 
form  of  the  unknown  pressure  loading  which  unfortunately  appears 

under  the  integral  sign.  As  was  done  in  the  subsonic  case  (Ref.  21), 
exact  theoretical  results  will  be  used  as  a weighting  function  and 
multiplied  by  an  assumed  functional  distribution  with  unknown  coefficients. 
The  form  of  the  weighting  function  is  particularly  important  in  supersonic 
flow  since  the  pressure  distribution  has  sharp  breaks  (discontinuities 
in  slope)  along  certain  Mach  lines.  This  procedure  has  been  used 
successfully  by  Cunningham  (Ref.  32).  The  form  of  the  pressure  loading 


used  herein  is 
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where 
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and 


~ I ® sin(m+l)0 
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in=0 


0 = Arcos(n) 


The  nondimens ional  chordwise  and  spanwise  variables  C.n  are 


(55) 


(56) 


C = 


c(ri) 


(57) 


and 


n = y/(b/2)  (58) 

The  weighting  function  P(x^,y^)  has  several  forms,  depending  on  whether 
the  wing  leading  and  trailing  edges  are  subsonic  or  supersonic.  The 
forms  used  will  be  based  on  exact  and  semi-exact  results  for  a general 
trapezoidal  wing  as  shown  in  Figure  15. 

For  a subsonic  leading  edge  and  supersonic  trailing  edge  (Fig  15a) , 
the  weighting  functions  used  are  the  pressure  equations  given  in 
Reference  40.  These  equations  are  based  on  an  exact  conical  flow  solution 
for  region  1 and  an  approximation  to  the  exact  results  in  the  tip  region  2. 
For  region  1: 


P(Xo,yo)  = 


(59) 


E(K)  o 
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where  E(K)  is  the  complete  elliptic  integral  of  the  second  kind  of 
modulus 

K = y[l  - b2/o2  (60) 

and  a = tan  A. 

In  region  2,  the  approximate  equation  from  Reference  40  is 
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(61) 


When  the  leading  edge  is  supersonic,  there  are  five  regions  to  be 
considered  (Figure  15b) . Exact  equations  for  each  region  are  given  in 
Reference  41,  which  is  based  on  Eward's  theory  (28). 

Region  1: 

P(x^,y^)  = = 4/-^s2-a2  (62) 

Region  2 : 


P(x 


.y^)  - [Arcos(f^)  + Arcos(§^); 


-OT 


(63) 
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Region  3 : 
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If  Che  trailing  edge  is  subsonic  (Figure  15c)  , the  pressure  functions 
calculated  above  for  either  case  (subsonic  or  supersonic  L.E.)  are 
corrected  by  a term  which  enforces  the  Kutta  condition  along  the  trailing 
edge.  This  term  is  an  approximation  based  on  the  exact  results  for  Region  3. 
It  is  used  in  the  same  manner  as  the  one  used  by  Cunningham  (Ref  32),  and 
is  multiplied  by  the  previously  calculated  pressure  terms: 


P(.Xo,yo)  = P(xE,y^) 
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where 
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The  function  P(x„,y  ) is  evaluated  at  the  y point  on  the  Mach  line  from 
E o o 

the  trailing  edge  vortex  using  the  appropriate  equation  from  (59)  through 

(68). 

With  the  desired  form  of  the  pressure  function  determined,  the  next 
task  is  to  evaluate  the  integral.  The  following  discussion  will  present 
the  necessary  equations.  The  evaluation  does  not  proceed  in  the  classical 
manner  of  treating  kernel  functions  (i.e.,  numerical  quadrature),  even 
though  functional  distributions  for  the  pressure  panels  from  the  discrete 
element  approach  (see,  for  example.  Woodward  (7))  is  used  to  simplify  the 
integration  procedure. 

First,  the  wing  surface  is  subdivided  into  a large  number  of  small 
finite  elements.  Equation  (53)  for  the  downwash  may  then  be  written 

w(x,y)  = I I Aw(x,y)  (74) 

S 

o 

where  Aw(x,y)  denotes  the  downwash  at  (x,y)  due  to  an  element  or  portion 
of  an  element  in  the  Mach  forecone  from  (x,y).  By  judiciously  choosing 
the  manner  in  which  the  wing  is  subdivided,  each  element  or  portion  there- 
of lying  in  the  region  of  integration  will  be  either  a rectangle  or 
triangle,  as  shown  in  Figure  16.  The  Integral  expression  for  a typical 
element  is 

'2 

AC  (x  ,y  ) (x-x  )dx  dy 

poo  o o o 

, : \2'  ' (7.5) 

^ ^o  y(x-x^)^-62(y-y^)^ 

Note  that  for  sufficiently  small  elements  the  pressure  loading  is 
essentially  constant  over  the  entire  element.  Consequently,  it  may  be 


Aw(x,y)  = 


^2 
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taken  outside  the  integral  and  evaluated  at  the  element  centroid  (x  ,y  ) 

o o 


with  the  result 


^2  ^^2 


4w(x,y)  = 
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The  integral  in  Equation  (76)  may  be  evaluated  in  closed  form  when  the 
appropriate  limits  of  Integration  are  known.  For  rectangular  elements 
such  as  element  A in  Figure  16 , the  result  is 


Aw(x,y)  = 


where 


K(Xo,yo)  = 


SAC  (x„,y  ) 


[-K(x2,y2)  + KCx^.y^^)  + K(xj^,y2)  -K(x^,yj^)]  (77) 
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For  triangular  elements  bounded  by  Mach  lines,  such  as  elements  B and  C, 

the  lower  x-limit  is  a function  of  y . The  form  of  this  function  causes 

o 

certain  terms  to  vanish  during  the  integration.  Results  for  both  elements 
B and  C are  then 


6AC  (x  ,y  ) 
Aw(x,y)  = E.— 2 — 2_ 

4lT 


[K(Xi,y2)  - K(x^,y^)] 


where  the  K-function  is  given  by  Equation  (78).  For  element  D,  which 
Includes  the  point  (x,y)  on  its  boundary,  two  separate  integrals  as  well 
as  the  poles  of  two  inverse-tangent  functions  must  be  evaluated.  The 
integration  and  evaluation  yields  the  well-known  Ackeret  result  (Ref.  34) 


for  an  infinite  span  wing: 


BAG  (x  ,y  ) 

Aw(x,y)  = poo 
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The  exact  integrals  for  elements  cut  by  the  wing  leading  edge,  such 
as  E and  F,  yield  three  complicated  expressions,  depending  on  whether  the 
leading  edge  is  subsonic,  sonic,  or  supersonic.  Rather  than  using  the 
complicated  results,  the  downwash  due  to  these  elements  is  approximated 
by  using  one-half  of  the  value  due  to  a rectangular  element  with  the  same 
X and  y dimensions.  The  outlines  of  such  elements  are  shown  by  dashed 
lines  in  Figure  16.  The  equations  for  the  downwash  is  then  one-half  the 
value  of  Equation  (77); 


BAG  (x  ,y  ) 
Aw(x,y)  = — 2_ 


[-K(x^y^)  +Kx2y^)  +K(x2y^)  +K(x^y2)  -K(x^y^) J (81) 


The  pressure  loading  evaluated  at  the  centroid  of  the 

triangular  element.  This  approximation  would  be  expected  to  have  little 
effect  except  for  control  points  close  to  the  leading  edge  of  very  highly 
swept  wings. 

The  only  remaining  task  is  to  evaluate  the  unknown  coefficients  B 

nm 

in  the  pressure  loading  Equation  (54) . This  requires  introducing  the 


tangency  condition 


sin  a(x,y)  + w(x,y)  = 0 


(82) 


where  sin  a(x,y)  is  the  freestream  velocity  component  normal  to  the  wing 
at  (x,y).  Equation  (82)  is  the  expression  for  no  flow  through  the  wing 
surface  and  must  be  satisfied  at  every  control  point  on  the  wing.  Using 
the  previously  developed  equations  leads  to 


H 


F 
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I I K^(x,y,x^,y^,X2.y2)  = - sin  a(x,y) 


where  is  the  appropriate  K-function  from  Equations  (77)  through 
(81)  (Note  that  = - tt  for  Equation  (80) . When  the  pressure  loading 
formula  Equation  (54)  is  inserted  into  Equation  (83),  there  are  (N+l)x(M+l) 
unknowns  in  the  resulting  eqpression.  Choosing  spanwise  and  chordwise 
control  points  according  to 

^m  = I , m=0,l,...,M  (8^) 


,n=0,l....,N 

j=n-hnN 


and  evaluating  Equation  (83)  at  each  control  point  will  give  a complete 

set  of  algebraic  equations,  one  for  each  unknown  coefficient  B . In- 

nm 

setting  the  pressure  formula  Equation  (54),  this  set  may  be  arranged  in 
the  matrix  form 


la^Kbj)  . 


- Y ^sin  a(x  ,y  ) } 
P j m 


where 


y sin{(nri-l)e} 


and  b,  is  the  column  vector  of  unknown  coefficients  B where  1 = n+mN . 
j nm 

Inversion  of  the  system  is  accomplished  with  a Gaussian  elimination 
algorithm. 
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Once  the  B coefficients  have  been  obtained,  the  value  of  the 
nm 

pressure  loading  is  knovm  at  every  point  on  the  wing  surface.  Values  of 
section  and  total  loads  and  corresponding  moments  may  then  be  obtained 
using  the  standard  aerodynamic  formulas  (Ref.  34)  and  numerical  quadrature. 

Results 

Before  attempting  to  calculate  loads  on  a warped  wing  in  the 
presence  of  a body,  it  was  felt  that  a numerical  study  was  necessary  to 
test  the  validity  and  accuracy  of  the  basic  method.  This  study  investi- 
gated the  effects  of:  (1)  number  of  and  placement  of  control  points; 

(2)  discretization  of  the  lifting  surface;  (3)  panel  size;  (4)  polynomial 
forms;  (5)  leading  edge  element  approximation;  and  (6)  flow  tangency. 

The  best  test  for  the  method  is,  of  course,  comparison  with  experiment. 
For  this  purpose,  two  low  aspect  ratio  wings  were  chosen,  one  delta  and 
one  trapezoidal.  Planform  schematics  and  (NC=6)x(NS=6)  control  point 
grids  for  each  wing  are  shown  in  Figure  17.  Theoretical  distributions  of 
pressure  loading  coefficient  AC^  were  compared  with  experimental  results 
(from  Reference  42)  at  6°  angle  of  attack  and  two  Mach  numbers,  1.61  and 
2.01. 

Some  representative  results  for  each  wing  are  shown  in  Figures  18 
through  21.  These  results  are  Intended  to  show  that  the  supersonic 
weighting  function  is  valid  and  that  the  general  method  of  solution  has 
good  agreement  with  theory,  even  for  the  highly  swept  delta  wing.  General 
results  and  some  of  the  problems  encountered  with  each  wing  will  be 
discussed  separately. 


The  first  example  in  Figures  18  and  19  is  the  aspect  ratio  1.342 
trapezoidal  wing.  The  wing  has  a 25°  leading  edge  sweep  and  a taper 


Figure  19.  Chordwlse  Pressure  Distribution  on  a Trapezoidal  Wing  in 
Supersonic  Flow  (M=«2.01). 
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Figure  21.  Chordwise  Pressure  Distribution  on  a Delta  Wing  in 
Supersonic  Flow  (M=2.01). 
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ratio  of  0.721.  The  airfoil  section  Is  a 3%  circular  arc,  w ch  Is 
essentially  a zero- thickness  wing  and  Ideally  suited  to  the  present 
theory.  The  leading  edge  Is  supersonic  for  both  Mach  numbers  shown. 

The  (NC*6)x(NS*6)  control  point  grid  shown  In  Figure  17  was  used  to 
obtain  the  results  shown.  The  wing  was  divided  Into  100  evenly  spaced 
spanwlse  panels  with  45  chordwlse  panels.  The  chordwise  panel  spacing 
used  the  cosine  distribution  which  placed  smaller  panels  near  the  leading 
and  trailing  edges.  Results  agree  exceptionally  well  with  experiment 
except  over  the  rear  half  of  the  chord  at  the  n=0.97  spanwlse  location. 
The  disagreement  In  this  region  Is  probably  due  to  viscous  effects. 

Figures  20  and  21  present  results  for  the  aspect  ratio  1.456  delta 
wing.  The  wing  has  a 70°  leading  edge  sweep,  which  Is  subsonic  for  both 
Mach  numbers  1.61  and  2.01.  The  chordwise  airfoil  section  is  a NACA 
65A003  thickness  distribution,  which  is  a 3%  thick  symmetric  airfoil  with 
a rounded  leading  edge.  The  (NC=6)x(NS”6)  control  point  grid  shown  in 
Figure  17  was  used  to  obtain  all  results  shown  for  this  wing.  Overall 
agreement  between  theory  and  experiment  la  quite  good,  the  major  dis- 
agreement occurring  in  the  expected  places  - near  the  blunt  leading  edge 
and  at  the  extreme  wing  tip. 

During  the  development  of  the  numerical  method  several  equations 
presented  in  the  theory  section  had  to  be  modified  to  obtain  acceptable 
results  for  the  delta  wing.  The  assumed  pressure  loading  form  used  for 
the  delta  wing  Is 
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(88) 


which  replaced  Equation  (54).  A cosine  distribution,  rather  than  even 
spacing,  was  used  for  the  spanwise  panel  size.  Smaller  panels  were 
placed  near  the  wing  tip.  The  control  point  placement  along  a chord  was 

i ^ (89) 

j=n+mN 

which  replaced  Equation  (85).  Numerical  studies  have  subsequently  shown 
that  these  changes  were  required  because  the  region  of  integration,  or 
forecone,  for  control  points  near  the  leading  edge  was  not  subdivided 
finely  enough.  This  subdivision  is  particularly  critical  for  highly 
swept  wings  with  subsonic  leading  edges. 

Numerical  tests  were  also  made  using  various  numbers  of  spanwise  and 
chordwise  panels  and  various  control  point  combinations.  Acceptable 
results  could  be  produced  with  a coarser  control  point  grid,  such  as 
(NC=3)x(NS=4) , and  in  some  cases  the  comparison  with  experiment  was 
slightly  better  than  that  shown  in  Figures  18  through  21.  The  (NC=6)x(NC=6) 
results  were  shown  in  the  current  report  for  two  reasons.  First,  it  is 
anticipated  that  for  a warped  wing  a large  number  of  control  points  will 
be  needed  to  accurately  define  the  pressure  field.  Second,  some  classical 
lifting  surface  methods  have  yielded  divergent  results  when  increasing 
numbers  of  control  points  were  used.  The  present  method  has  not  exhibited 
such  trends,  even  for  (NC“6)x(NS“8)  grids.  Typical  computation  times  were 
45  seconds  for  the  (3x4)  case  and  118  seconds  for  the  (6x6)  case  on  an 
IBM  370/155. 


Conclusions 

The  original  intent  during  development  of  the  supersonic  program  was 


to  use  the  existing  subsonic  program  and  make  the  two  almost  Identical  - 
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the  only  mofifications  being  the  closed  form  Kernel  integrals  and  the 
limiting  region  of  integration  for  a control  point  from  the  entire  wing 
to  the  forecone.  It  was  felt  that  this  approach  would  require  minimal 
effort  and  allow  rapid  development  of  the  supersonic  method.  Such  was 
not  to  be  the  case,  however. 

It  became  immediately  apparent  that  different  weighting  functions 
were  required  for  the  assumed  pressure  loading,  and  considerable  time 
was  required  to  develop  these  from  the  exact  forms  given  in  the  liter- 
ature (Refs.  40,  41,  43).  However,  as  Cunningham  discovered  (Ref.  32), 
results  using  these  functions  were  far  superior  to  those  using  the  sub- 
sonic program  forms. 

Several  other  conclusions  which  need  to  be  noted  were  reached  during 
the  numerical  study  of  the  basic  supersonic  method.  In  the  current 
version,  the  placement  of  control  points  must  vary  depending  on  whether 
the  leading  edge  is  subsonic  or  supersonic.  The  area  of  integration  for 
each  control  point  needs  to  be  subdivided  better,  particularly  for  con- 
trol points  near  the  leading  edge.  The  latter  change  may  allow  standard- 
ization of  control  point  placement,  and  may  also  minimize  the  effects  of 
the  leading  edge  element  approximation,  even  on  highly  swept-wlngs.  A 
new  form  for  the  assumed  polynomial  distribution  is  indicated.  This  new 
form  should  be  applicable  to  wings  of  any  taper  ratio.  This  need  may 
become  even  more  apparent  when  body  effects  and  wing  warpage  are  Included. 

Finally,  planar  wing  results  using  the  method  in  its  current  form 
show  good  agreement  with  experiment  and  are  not  divergent  with  increasing 
numbers  of  control  points.  It  is  felt,  however,  that  any  or  all  of  the 
previously  mentioned  changes  could  noticeably  Improve  the  results  and 
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possibly  the  convergence  of  the  method.  It  would  also  be  desirable  if 
results  and  computation  requirements  could  be  compared  with  other  numer- 
ical methods  such  as  References  31,  7,  or  32. 
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APPENDIX 

DERIVATION  OF  THE  COMPRESSIBLE  LIFTING  SURFACE  POTENTIAL 
The  following  derivation  parallels  that  given  in  Reference  3 for 
incompressible  flow,  which  begins  with  the  most  elementary  solution  of 
the  linearized  potential  equation,  the  point  source.  The  potential  for 
a point  source  in  compressible  flow  is 


(A-1) 


where  m is  the  source  strength,  r is  the  elliptic  radius  defined  by 


r = 


and  the  Mach  parameter  6 is  /l-M2  . 

The  potential  for  a doublet  with  its  axis  oriented  in  the  z-direction 
is  found  by  differentiating  (A-1)  with  respect  to  z: 
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The  solution  for  the  elemental  horseshoe  vortex  is  the  integral  of  a 
semi-infinite  line  of  doublets: 
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Performing  the  integration  and  evaluating  the  limits  gives: 
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which  differs  from  the  Incompressible  results  only  by  the  6^  terms  in  the 
radical. 

Now  consider  a sheet  of  these  elemental  vortices,  lying  in  the  z=0 
plane.  The  constant  m is  replaced  by  Y(*Q»yQ)  the  vorticlty  per  unit 
area,  which  is  assumed  to  be  a function  of  the  position  on  the  surface. 
The  potential  at  any  arbitrary  point  (x,y,z)  in  the  flow  field  due  to  an 
infinitesimal  area  dx  dy  of  the  sheet  is 
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dx  dy  (A-6) 
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Denoting  the  total  surface  area  of  the  sheet  by  S,  the  potential  due 
to  the  entire  surface  is  then 


i}i(x,y,z)  = ffdHx  ty.z) 
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For  supersonic  flow,  the  result  is 
'>’(x,y,z) 


i-  f f 

2^  J J fv-i,  ■»  24.-2 


z (x-x  ) 
o 


^o^  V(x-x^)^-e^ (y-y^)^-6^z2 


dx  dy  (A-9) 
o o 


